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Abstract 



We present measurements of the production of Ds~ mesons relative 
to Dg'^ mesons as functions of Xi? and of pt^ for a sample of 2445 Dg de- 
cays to (fm. The Dg mesons were produced in Fermilab experiment E791 
with 500 GeV/c Tr~ mesons incident on one platinum and four carbon 
foil targets. The acceptance-corrected integrated asymmetry in the xp 
range -0.1 to 0.5 for D^^ mesons is 0.032 ± 0.022 ± 0.022, consistent 
with no net asymmetry. We compare the results as functions of Xi;' and 
Pt^ to predictions and to the large production asymmetry observed for 

mesons in the same experiment. These comparisons support the 
hypothesis that production asymmetries come from the fragmentation 
process and not from the charm quark production itself. 
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Previous studies of charm particle production in hadron beams have shown 
large enhancements in the forward production of charm particles that contain a quark 
or di-quark in common with the beam (leading particles) relative to those that do not 
(non- leading particles). Neither leading-order nor next-to-leading-order perturbative 
QCD calculations of charm quark production can account for the observed asymme- 
tries |[TO| . |TT| . The phenomenon is due either to unexpected contributions to charm 
quark production or to features of the hadronization process. 

Three classes of models have been proposed to account for hadronization in the 
beam fragmentation region: coalescence of produced charm quarks with valence 
quarks or diquarks from the beam [p!^ -[l^, coalescence of charm and valence quarks 
when both originate in the beam particle 0-0, and string fragmentation as im- 



plemented in the Lund Model of the PYTHIA Monte Carlo program |118[ and by 
Piskounova |TP|. All these models involve a mechanism of attachment between the 
produced charm quarks and the valence quarks in the beam, and include a varying 
amount of actual coalescence. This coalescence is the dominant mechanism for the 
leading particle effect in the first two classes of models. They include no mechanism 
which will lead to a significant production asymmetry between particles except when 
one of them has a quark in common with the beam and the other does not. Since 
the probability of coalescence between two quarks is larger when they have similiar 
velocities, the production asymmetry between leading and nonleading charm mesons 
is expected to grow larger with increasing xp and smaller with increasing pt^ . While 
data from Fermilab experiment E791 do show the expected increase in asymmetry 
with increasing xp-, they do not support the prediction of decreasing asymmetry with 
increasing 0. 

In string fragmentation models, specifically the PYTHIA/LUND implementation, 
a produced charm quark is always initially connected via color strings to valence 
quarks or diquarks in the beam and the target. A leading particle can be produced 
by a meson beam in cases where the string attaches a charm quark to a beam anti- 
quark and/or an anticharm quark to a beam quark. Some fraction of the time, the 



3 



string invariant mass is too small to allow the production of a quark- ant iquark pair, 
and the process reduces to that of final state coalescence, with the same general 
features. However, in most cases additional qq pairs are produced in the evolution of 
the strings. The string fragmentation model predicts an asymmetry even at relatively 
small values of xp- This asymmetry is predicted to be relatively flat as a function 
of j9t^, consistent with E791 measurements [0. String fragmentation also allows for 
a possible asymmetry in cases where there is no quark in common with the beam, if 
one of the quarks is produced in the string fragmentation. 

For the nearly pure n~{ud) beam pO[ in Fermilab experiment E791 ||23|, D^{cd) is 
a leading particle and D^{cd) is the corresponding nonleading particle. We have pre- 
viously reported a significant production asymmetry, increasing with xp and rel- 
atively flat in pt^, and well represented by a modified version of the PYTHIA/LUND 
Monte Carlo [0]. In contrast, Z)+(cs) and D~{cs) mesons do not share any valence 
quark type with the pion beam. There should thus be no asymmetry between Ds~ 
and production if the effect is due to coalescence. The limited Dg sample sizes in 
previous experiments did not allow a detailed measurement of the production asym- 
metry for Dg mesons PJ2^. Indeed, no asjTiimetry data has been reported for 
pion-induced production. We report the first results on the Dg production asymme- 
try in pion-nucleon collisions, and with sufficient data to study the asymmetry as 
functions of xp and of pt^. 

The E791 spectrometer had 23 silicon microstrip detector planes (SMD's); six 
upstream of the target foils, seventeen downstream, to measure charged particle tra- 
jectories. This SMD system allows precise measurement of the position of the charm 
production vertex (primary vertex) and charm decay vertex (secondary vertex). Two 
magnets, together with the downstream SMD's, 35 drift chamber planes, and two 
proportional wire chamber planes, provided momentum measurement. Particle iden- 
tification of pions, kaons, and protons was obtained with two segmented, threshold 
Cerenkov counters. The total E791 recorded data sample is 2 x 10^'' 500 GeV/c tt~ 
nucleon interactions, which satisfied a loose requirement on the total transverse en- 
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ergy deposited in the calorimeter, assuming that particles originate at the target foils. 
After data reconstruction, only those interactions in which one or more decay vertex 
candidates were detected were kept for further analysis. 

From the above sample, charm meson decays to (f)7i were sought by starting with 
candidate 4> K^K~ decays. The (j) selection began with combinations of two 
oppositely charged, well-reconstructed tracks, at least one of which was identified as 
a kaon by the Cerenkov counters. The combination was retained as a candidate 
if the invariant mass of the two tracks, assumed to be kaons, was within 10 MeV/c^ 
of the (f) mass (1.0194 GeV/c^). For each candidate, a third track, assumed to be 
a pion, was sought, for which the invariant mass of the three tracks, assumed to be 
KKir, was within a mass window 1.79 - 2.05 GeV/c^, which includes both the Ds 
and the D"^ masses. A fit to the point of intersection of the three tracks was then 
performed. 

To further distinguish charm meson decays from background, additional selection 
criteria were applied to each secondary vertex candidate. The values of the selection 
criteria were chosen to maximize the statistical significance {Ns/ ^(TvT-K/Vb)) for each 
xp bin, where Ng is the number of Dg signal events and is the number of back- 
ground events under the Ds mass peak. The signal is represented by fully simulated 
Dg events generated by the PYTHIA Monte Carlo and the background is E791 data 
for which the candidate mass is in the sideband of the Dg peak, but outside the 
signal region. For most of the xp range, typical cuts are as follows. The longitudinal 
separation between the primary and secondary vertices was required to be more than 
10 times the experimental resolution for the measured separation. The impact pa- 
rameter of the reconstructed momentum vector of the Dg candidate with respect to 
the primary vertex was required to be less than 40 microns. The ratio of the distances 
of closest approach to the secondary and primary vertices were computed for each 
of the three decay tracks. The product of the three ratios was required to be less 
than 0.006. This means that in most cases the tracks that make the secondary vertex 
should be significantly closer to the secondary vertex than to the primary vertex. 
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It was required that no track in the event, except the three making the secondary 
vertex, pass within 10 microns of the secondary vertex. 

Figure 1 shows the K^K^t^^ and K^K^tt^ invariant mass distributions after 
the above selection criteria are apphed. A fit with two Gaussian curves and a hnear 
background apphed to the combined spectrum yields 2445 ± 60 D^^ mesons and 1449 
± 55 mesons. Fits to the individual spectra yield 1188 ± 41 D^^niesons, 1257 ± 
42 Ds~ mesons, 638 ± 37 mesons and 811 ± 40 D~ mesons. These uncorrected 
data show a clear excess of D~ mesons over mesons, consistent with the earlier 
measurement by this collaboration for the more copious decay D"^ K^-k^-k^ 0, but 
show no statistically significant asymmetry in the numbers of Ds^ and Dg^ mesons. 

In order to study the asymmetry in the production of Dg^ mesons, K^K^tx^ 
invariant mass plots were made for particular ranges oi xp and pt^. The yields of Dg^ 
and Ds~ mesons were then measured by fitting the peaks in the mass plots to Gaussian 
curves of fixed width and mean plus a linear background, and then integrating the 
area under the Gaussian curves. The evolution of the width of the Gaussian curves 
with xp and pt^ was determined from Monte Carlo simulation, and is the same for 
both charge states. 

An asymmetry parameter A for Dg production is defined as 

AV^AT^ 

Here, A^^^- is the number of negatively charged Dg mesons and Nj~i^+ is the number 
of positively charged Dg mesons produced in a particular bin of or pt^ . 

From studies using our complete Monte Carlo simulation, we find that our accep- 
tance, which includes the geometrical acceptance as well as the reconstruction effi- 
ciency, falls off rapidly at high xp- This inefficiency at high xp increased with time 
during the data taking period. This was mainly due to beam effects which caused a 
cumulative local loss of efficiency in the drift chambers. A correction ranging from 
—6% to +12% was applied bin-by-bin for the difference between the acceptances of 
positive and negative Dg mesons. Residual uncertainties in this correction are treated 
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as systematic errors and are taken from reference 0, using the r.m.s. spread of four 
Monte Carlo simulations representing different time periods during the data taking. 

A second source of systematic uncertainty is due to mesons produced by the 
2.5% K~{us) content of the incident beam for which the integrated asymmetry A 
is measured by E769, at 250 GeV/c beam energy, to be 0.25 ±0.11 No correction 
to the asymmetry is applied, but we assign a systematic error of 0.006 in the xp 
range (—0.1 to 0.3) and 0.015 in the xp range (0.3-0.5) to account for this source of 
systematic uncertainty. The increase in systematic error at larger xp accounts for 
greater uncertainty in the K~ contribution in the forward direction. 

The acceptance-corrected production asymmetry for Ds mesons in the xp range 
—0.1 to 0.5 is 0.032 ± 0.022 ± 0.022, consistent with no net asymmetry. Here the 
first error is the statistical error and the second error is the systematic error from the 
sources described above. 

The acceptance-corrected asymmetries as functions oi xp and pt^ are shown in 
Tables I and II and in Figures 2 and 3 respectively. A line at zero asymmetry is 
drawn in the two figures for reference. Asymmetries in the production of Ds~ and 



Dg^ mesons as predicted by a tuned version of the PYTHIA/LUND Monte Carlo ^1 



are also shown in Figures 2 and 3 for comparison to the measured asymmetries. In 
this version of PYTHIA, parameters such as the mass of the charm quark and the 
intrinsic transverse momentum of the partons have been tuned so that the asymmetry 
of the Monte Carlo sample matches the asymmetry for the — ^ Kt^tx using E791 
data 0]. While E791 data are consistent with the PYTHIA model in both xp and 
j9t^, a better is obtained for the hypothesis of A=Q in both cases. Figures 2 and 
3 also show a comparison of Ds and asymmetries as a function oi xp and pt^ 
measured using E791 data. The data are from a high statistics measurement 
using the decay mode Kn-n The production asymmetry is larger for D"^ 

mesons than for Ds mesons for all Xp and pt^ . An intrinsic charm model predicts zero 
asymmetry for Ds production at all x p and pt^ |jl6[ . 



In summary, we have observed both and Df decays in the same apparatus and 
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find a net production asymmetry only for the D . We find that the Dg asymmetry 
is consistent with zero in each bin oi xp and pt^ in the ranges covered by E791. Since 
the beam used in experiment E791 contains almost no strange particles and the Ds 
asymmetry is consistent with zero, it can be concluded that the asymmetry seen in 
the D"^ mesons is not likely to come from any asymmetry in the production of the 
charm quarks themselves, but originates in the hadronization process. 
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is operated by the Universities Research Association, Inc., under contract with the 
United States Department of Energy. 



8 



REFERENCES 

[1] E595 Collaboration, J.L. Ritchie et al, Phys. Lett. B138, 213 (1984). 

[2] NA27 Collaboration, M. Aguilar-Benitez et a/., Phys. Lett. B161, 400 (1985). 

[3] NA32 Collaboration, S. Barlag et al, Z. Phys. C49, 555 (1991). 

[4] WA82 Collaboration, M. Adamovich et al, Phys. Lett. B305, 402 (1993). 

[5] E769 Collaboration, G.A. Alves et al, Phys. Rev. Lett. 72, 812 (1994). 

[6] E769 Collaboration, G.A. Alves et al, Phys. Rev. Lett. 77, 2388 (1996). 

[7] E791 Collaboration, E.M. Aitala et al, Phys. Lett. B371, 157 (1996). 

[8] WA89 Collaboration, R. Werding, Proceedings XXVII Int. Conf. on High Energy 
Physics, Glasgow, U.K., 20-27 July, 1994, p. 1023. 

[9] Beatrice Collaboration, M. Adamovich et al, CERN Preprint CERN-PPE/96- 
180 (1996). 

[10] P. Nason, S. Dawson and R.K. Ellis, Nucl. Phys. B327, 49 (1989). 

[11] S. Prixione et al, Nucl. Phys. B431, 453 (1994). 

[12] B.L. Combridge, Nucl. Phys. B151, 429 (1979). 

[13] V. Barger, F. Halzen, and W.Y. Keung, Phys. Rev. D25, 112 (1982). 

[14] P. Mazzanti and S. Wada, Phys. Rev. D26, 602 (1982). 

[15] S.J. Brodsky et al, Phys. Lett. 93B, 451 (1980); R. Vogt and S.J. Brodsky, 
Nucl. Phys. B 438 261 (1995). 

[16] R. Vogt and S.J. Brodsky Nucl. Phys. B478, 311 (1996). 

[17] R.C. Hwa, Phys. Rev. D51, 85 (1995). 

[18] H.-U. Bengtsson and T. Sjostrand, Comp. Phys. Comm. 46, 43 (1987). 

9 



[19] 0.1. Piskounova, Nucl. Phys Proc. Suppl. 50, 179 (1996). 
[20] A. Malensek, private communication. 

[21] Specific parameters changed in the PYTHIA Monte Carlo (version 5.7 with ver- 
sion 7.3 of JETSET) were PMAS(4.1) (the c quark mass in GeV/c^) changed 
from 1.35 to 1.7, PARP(91) (the average kt'^ in (GeV/c)^) increased from 0.44 to 
1.0, PARP(93) (the maximum allowable kt in (GeV/c)) increased from 2.0 to 5.0, 
and MSTP(92) (the remnant quark-diquark energy sphtting function) changed 
from 4 to 3. 

[22] E769 Collaboration, G.A. Alves et ai, Phys. Rev. Lett. 77, 2392 (1996). 

[23] J.A. Appel, Ann. Rev. Nucl. Part. Sci. 42, 367 (1992) and references therein; 
D.J. Summers et al, Proceedings of the XXVII^'^ Rencontre de Moriond, Elec- 

trowcak Interactions and Unified Theories, Lcs Arcs, France (15-22) March 1992, 
417; S. Amato et a/., Nucl. Instr. Methods A, 324, 535 (1993). 



10 



TABLE I. Dg production asymmetry as a function oi xp from E791 data, integrated 
over the pt' range to 10 [GeV/cY. The first error is statistical and the second 
systematic. 





Dg asymmetry A 


-0.1 


- 0.0 


0.015 ± 0.054 ± 0.023 


0.0 - 


- 0.05 


-0.002 ±0.045 ±0.007 


0.05 


- 0.1 


0.001 ± 0.047 ± 0.020 


0.1 


- 0.2 


0.020 ± 0.042 ± 0.027 


0.2 - 


- 0.25 


0.080 ± 0.091 ± 0.018 


0.25 


- 0.3 


0.083 ± 0.125 ± 0.036 


0.3 


- 0.4 


0.104 ± 0.140 ± 0.075 


0.4 


- 0.5 


0.149 ± 0.252 ± 0.105 



TABLE II. Dg production asymmetry as a function of pt^ from E791 data, inte- 
grated over the xp range —0.1 to 0.5. The first error is statistical and the second 
systematic. 



Pt\GeY/cf 


Dg asymmetry A 


0.0 - 1.0 


-0.016 ±0.035 ±0.013 


1.0 - 2.0 


0.036 ± 0.043 ± 0.036 


2.0 - 3.0 


0.089 ± 0.052 ± 0.056 


3.0 - 4.0 


0.081 ± 0.056 ± 0.089 


4.0 - 7.0 


0.032 ± 0.053 ± 0.078 


7.0 - 10.0 


-0.077 ± 0.053 ± 0.068 
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FIG. 1. K^K TT and K^K tt"*" invariant mass plots. 
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FIG. 2. Comparison of the Dg (solid circles) production asymmetry as a function of 

Xi? to the (open circles) production asymmetry, both measured by experiment E791. 
The acceptance-corrected data are integrated over the pt^ interval (1-10) (GeV/c)^. Also 
shown in the figure is the Dg production asymmetry (solid line) predicted by the "tuned" 
PYTHIA Monte Carlo described in the text. 
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FIG. 3. Comparison of the Dg (solid circles) production asymmetry as a function of pt'^ 
to the (open circles) production asymmetry, both measured by experiment E791. The 
acceptance-corrected data are integrated over the xp interval (—0.1, 0.5) for Dg mesons and 
(-0.1, 0.8) for D"^ mesons. Also shown in the figure is the Dg production asymmetry (solid 
line) predicted by the "tuned" PYTHIA Monte Carlo described in the text. 
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